MicroChannel flow of a shear-banding fluid: enhanced conflnement effect and 

interfacial instabihty 
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Using a micro particle imaging velocity technique, we resolve for the first time the three dimen- 
sionnal structure of wormlike shear banding flows in straight microchannels. The study revealed 
two effects, which should be generic for shear banding flows: the first is a strong amplification of 
the confinement induced by the edge of the channel, the second is an instability of the interface 
between the shear bands. A detailed quantitative comparison of our experimental measurements 
with a theoretical study of the diffusive Johnson Segalman model leads to excellent agreement. Our 
study clarifies the nature of shear banding flow instabilities, and shows that, despite the challenging 
complexity of the situation and the uncertainty regarding their molecular structure, shear banding 
flows in confined geometries are amenable to quantitative modelling, a feature that opens pathways 
to their practical utilization. 



In many complex fluids, a strongly nonlinear coupling 
between flow and microstructure can trigger the onset 
of flow instabilities, leading to the formation of coexist- 
ing "shear bands" that support differing shear rates at 
a common shear stress. This model situation for out- 
of-equilibrium phase separation has been the focus of 
intense research in recent years with wormlike mi- 
cellar surfactants the most widely studied systems [Ij, 
both as model materials and for their widespread appli- 
cation in, e.g.^ oil recovery and personal care products. 
Significant progress has been achieved by imaging experi- 
mentally the local flow profile, first by NMR [3] and later 
by DLS [4], ultrasound or particle image velocime- 
try [6|. In this Letter we report a significant advance 
in such efforts, with the first full three dimensional (3D) 
velocity imaging of a shear-banding fluid. We perform 
this in the context of pressure driven flow in a rectilinear 
microchannel of rectangular cross section. 




Using this imaging, we report two flow phenomena. 
The first is a strongly enhanced confinement effect in- 
duced by the edges, that we believe may prove generic in 
a broad class of heterogeneous confined flows. The second 
is an instability of the interface between the shear bands, 
characterized by a downstream evolution from (near the 
inlet) a smoothly bowed interface between the bands, to 
(far downstream) a modulated interface with wavevector 
in the vorticity direction. Throughout, we compare our 
data with a theoretical study of the diffusive Johnson- 
Segalman (DJS) model in the same geometry. The ex- 
cellent quantitative agreement we find without fitting 
assesses the robustness of our observations in the gen- 
eral context of shear-banded flows. It also suggests that 
the behavior of wormlike micelles in complex geometries 
is amenable to theoretical description at a quantitative 
level. 

In particular, the interfacial instability we observe re- 
lates to a growing body of experiments showing com- 



FIG. 1: Squares: experimental steady state flow curve [ill ]. 
Line: theoretical curve with the DJS model for a — 0.3, ?7 = 
0.05. Both are renormalized by the shear rate 7c at the onset 
of shear banding, and by the value of the stress plateau Gc- 
Correspondence is not expected at high shear rates (see main). 
Inset: sketch of the basic (ID) shear-banded flow, with the 
coordinate system used throughout the paper. 



plex spatio-temporal patterns in shear-banded flow, most 
of which have been to date in curved Couette geome- 
try. While observations of instabilities are now quite 
widespread in that geometry via ID measurements, 
the only previous explicit observation of a 2D modula- 
tion of the interface is in Ref.[7]. The theoretical un- 
derstanding of this instability remains preliminary 
but suggestions about its origin include (i) a bulk vis- 
coelastic Taylor- Couette-like instability originating in the 
strongly sheared band and (ii) an interfacial insta- 
bility triggered by the normal stress jump between the 
bands [9|, llO|. An important additional contribution of 
this Letter is to give evidence favouring (ii), by eliminat- 
ing the cell curvature needed for (i). 

Experimental — The sample comprises CTAB 0.3M, 



NaNOs 0.405M (Acros organics), which forms a semi- 
dilute solution of wormlike micelles. Its basic rheol- 
ogy was characterized in Ref. [11], with a stress plateau 
of 105 Pa beginning at 7c = 5 corresponding to 
a pressure gradient for the onset of shear-banding of 
Gb = 3.10^ Pa. m~^. Our experimental setup com- 
prises pressure driven flow in a microchannel of rect- 
angular cross section, with dimensions {Lx^Ly^Lz) = 
(5 cm, 64/xm, 1 mm) with the axes of Fig. [H The pho- 
tocurable glue used in the fabrication permits high aspect 
ratio (1 : 16) channels that do not deform even under 
pressures far in excess of those in this study • In this 
work, the setup of [11] has been substantially modified 
to allow full 3D velocity imaging, as follows. Fluorescent 
tracers of diameter 500 nm are tracked at 50 Hz with a 
charged coupled device camera (Allied) through a lOOx 
oil immersion objective mounted on a piezo, giving a 1/im 
thick focal plane that can be displaced in y with 10 nm 
resolution. Images comprise (Ax, Az) = {SG/j^m, 72/im) 
windows, the position of which is determined using a 
moving table (Marzhauser). These images can be divided 
into subimages along z, resulting in correlations between 
images of size 9 /im for z. As a whole, then, velocity is 
mapped over the sample with spatio-temporal resolution 
(Ax, Ay, Az, At) = (72 /im, 1 /im, 9 /im, 40 ms). 

Theoretical — We compare all our data with the pre- 
dictions of the diffusive Johnson Segalman (DJS) model 
. This is not intended to give a microscopically faith- 
ful description of our particular fiuid, but is instead the 
minimal phenomenological tensorial model to capture the 
constitutive non-monotonicity needed for banding. Its 
parameters are: solvent viscosity r]; slip parameter a; lin- 
ear modulus Go; relaxation time r; and interfacial thick- 
ness normalized by the height l/Ly. Values for these 
are guided by previous linear and nonlinear bulk rheo- 
logical measurements on this system, independent of the 
confinement effect or instability studied here. Specifi- 
cally, the non-linear fiow curve of Fig. [T] suggests a ratio 
Tj/G^T ~ 0.05 of high and low shear phase viscosities. 
Perfect correspondence is impossible because the DJS 
model is oversimplified in having a Newtonian high shear 
branch, compared to the sublinear experimental trend. 
For pressure drops quite close to the onset of banding, as 
here, we do not expect this limitation to be severe. We set 
a = 0.3, though our results are robust to variations in this 
quantity. We estimate / ^ (^bT/Gq^^)^/^ = O(lO-^m), 
based on linear rheology. The experimental aspect ra- 
tio is Lz/Ly = 16; numerically we explore 8,16,32,64. 
Numerical calculations are peformed initially in units in 
which Go = l,T = l, L^ = l, but results are usu- 
ally presented by scaling stresses (or pressure drops) and 
shear rates by their values at the onset of banding, as in 
Fig. 1. Below we will demonstrate excellent quantitative 
agreement between theory and experiment for both the 
confinement effect and the interfacial instability, without 
any fitting specific to these phenomena. 
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FIG. 2: (a) Experimental velocity profiles across the centre- 
line y — Ly/2 of the cell cross section in the zone near the 
inlet (a typical error bar of 5% is not shown for clarity). Pres- 
sure drop G /Gh — 0.9, 1.06, 1.3, 1.6 (solid, dotted, dashed and 
dot-dashed lines) with Gh = 3.10^ Pa.m~^; (b) correspond- 
ing numerical results for the same scaled pressure drops with 
a = 0.3, T] = 0.05, / = 0.01. (c) Experimental (circles) and 
theoretical (squares) penetration length versus scaled pres- 
sure drop, (d) Steady state greyscale snapshot of invariant 
shear rate \\/vx \ for pressure drop G/Gh = 1.02, 1.09, 1.30 for 
Lz/Ly = 8.0 (aspect ratio rescaled for clarity). 



Our numerical study considers flow driven along a 
channel by a constant pressure drop G = —dxP-, assum- 
ing translational invariance of the velocity field in the 
main flow direction x. It comprises two separate parts, 
to allow a cleaner understanding of the two flow phenom- 
ena seen experimentally. First we study the enhanced 
conflnement effect by simulating channels of high aspect 
ratio with closed walls in y and z. In order to isolate 
this conflnement effect from the interfacial instability, in 
this part of the study we disallow any secondary flows in 
the channel cross section, thus artiflcially switching off 
the instability. Our aim is thereby to capture the "basic 
state" seen experimentally close to the inlet, before the 
instability sets in downstream. The second part of our 
study focuses on the interfacial instability, which we now 
isolate from the conflnement effect by switching off the 
latter via the use of periodic boundaries in z. This work 
is thereby intended to approximate the flow seen exper- 
imentally in the central region of z, well away from the 
lateral walls. In each run we evolve the code to steady 
state t ^ oo with the aim of capturing the flow pattern 
far downstream from the inlet x ^ oo, assuming that 
these limits correspond to each other. A stringent check 
of our code, adapted from a former study [M], is provided 
by comparing against a true linear stability analysis its 
linear regime growth dynamics of tiny amplitude pertur- 
bations to an initially flat interface. 
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FIG. 3: Greyscale of the experimentally measured velocity 
component Vx in the channel cross-section for G/Gh = 1-3, at 
an early downstream location x = 1 cm before the instability 
has set in (a) and far downstream x = 4 cm (b), with the 
fully developed pattern, (c) Corresponding numerical snap- 
shot of the fully developed pattern with G/Gh = 1-3, a = 0.3, 
Tj = 0.05, / = 0.0015, Lz = 6.0 and periodic boundaries in z. 
(d) Numerically computed secondary flow field for the same 
conditions. 



Results I: basic state and confinement effect — First 
we explore the "basic state" that is observed across the 
channel cross-section in the zone far enough from the in- 
let for the local rheology to have attained a stationary 
banded state, but close enough that the interfacial insta- 
bility has not set in, paying close attention to the role of 
the lateral walls of the channel {z = {0, Lz}) in break- 
ing translational invariance in the z direction. Given the 
high aspect ratio of our channel {L^/Ly = 16), we might 
expect the flow profile to be almost z independant except 
for regions close to the lateral walls. Indeed for Newto- 
nian flows it is known that for ^ Ly these regions 
have a size of order Ly. The new phenomenon we re- 
port is that for pressure gradients G beyond the onset 
of shear-banding, the velocity field can vary significantly 
with z across the entire channel cross-section. 

To quantify this, we measure for each G the cen- 
terline velocity profiles Vc{z) = Vx{y = Ly/2^z) and 
a penetration length A;^, characterising the size of the 
lateral regions affected by the walls. is defined by 
Vc{^z) = 0.96 with Vm = Vc{z = Lz/2) the maximum 



velocity. This criterion is chosen such that A^ = Ly for 
Newtonian flow for Lz ^ Ly. For our shear banding 
fluid, Vc{z) profiles (Fig. [2^,b) and A^ plotted vs. pres- 
sure gradient (Fig.[2t) compare with striking quantitative 
agreement between theory and experiment. 

Three regimes are evident. Below the onset of band- 
ing (regime I, G/G\y < 1) the flow is Newtonian to good 
approximation: Vc{z) profiles are invariant over a large 
zone around the centrepoints Lz/2 (Fig. [2^,b solid line) 
and Kz ^ Ly as remarked above. At the onset G = Gb of 
regime II a thin band of high shear rate nucleates along 
each of the long walls y = {0, L^} around the centre- 
points z = Lz/2, these being the locations at which the 
shear stress is maximum and so first exceeds the plateau 
(banding) stress of the bulk flow curve as G is increased. 
These bands are seen in our numerical y — z maps of the 
invariant shear rate \ Vvx \ in Fig.[2]i, top. It is at this on- 
set that the centreline profile suddenly becomes strongly 
bowed (Fig.[2^,b, dotted line). Reflecting this is the dra- 
matic jump of Az to Az ^ 5 Ly. As G increases further 
the thickness of the high shear bands increases, and so 
does their lateral spread along the long walls (Fig. [2]i, 
middle). At the onset of regime III, shear bands finally 
form along the short walls (Fig.[2]l, bottom; and Fig.[2^,b 
long dashed and dot-dashed). Az then shows a downward 
kink into a decline to Az ~ Ly at high G. 

How can the large value of Az in regime II be under- 
stood? For a shear-banding fluid, smooth stress inho- 
mogeneities across the channel cross section are strongly 
magnified in the velocity field via the nucleation of lu- 
bricating bands at the long walls. For a complete slip at 
these walls the system would "forget" the short dimen- 
sion altogether and show parabolic Poiseuille flow across 
the long centreline. Simulations at fixed Lz in regime II 
indeed show Az to increase with decreasing viscosity of 
the high shear band (not shown). This has important im- 
plications for a broad class of flows with stress dependent 
wall slip, to be explored in future. 

Results U: instability — ^We describe first the steady 
state experimental velocity profiles Vx{y^z) at different 
downstream locations x, for a fixed pressure drop, in a 
cross sectional area [z = — 1 mm] by [?/ = — 14 /im] 
containing the high shear phase together with the inter- 
face near the bottom wall. Close to the inlet {x = 1 cm, 
Fig. [3^) the interface between the bands is smoothly 
bowed, corresponding to the basic state discussed above. 
In contrast, far downstream (x = 4 cm. Fig. [Sb) the in- 
terface has undulations with wave vector along the z axis, 
with wavelength of the order of the channel height Ly. 
This pattern is steady in time, and invariant with re- 
spect to further progression downstream. (At intermedi- 
ate downstream locations 1cm < x < Acm, the instability 
develops via interfacial corrugations that emerge close to 
the lateral walls, and progressively invade towards the 
cell centre with increasing x; not shown.) 

Such interface instabilities have already been predicted 
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FIG. 4: Solid circles: experimental wavevector of the interfa- 
cial undulations as a function of applied pressure drop. Open 
squares: corresponding numerical results for the fully non- 
linear state computed with 77 = 0.05, / = 0.0015, Lz — 4.0. 
The numerical wavevector is constrained by the cell length to 
be heavily quantised as qz = 2m^/Lz. We therefore also show 
results (open circles) for the maximally unstable mode in a 
linear stability analysis, which does not suffer quantisation. 
Insert shows the experimental amplitude A in units of Ly. 



[iqI , [i4 | and observed developing in time in Couette 
geometry. Similarly to these former predictions, ongo- 
ing linear stability analysis of Poiseuille flow shows the 
concurrence of unstable Qx and Qz interfacial modes of 
wavelengths comparable to Ly , over a range of parameter 
values. Unfortunately our measurement averages veloci- 
ties over Ax = 72 jim > Ly which precludes observation 
of wavevectors Qx. We do, though, concentrate on explor- 
ing with our setup any instability to modulation in the 
vorticity direction. Thanks to our flow cell aspect ratio 
Lz/Ly = 16, there is a domain of z around the middle 
of the channel {z = Lz/2) over which the curvature of 
the basic state due to the enhanced confinement effect 
(Fig. [3tt) is negligible compared to the interface modula- 
tion for around ten wavelengths (Fig. [8)3). We perform 
our analysis in this large aspect ratio domain that can 
be considered invariant in z. Our corresponding numer- 
ics are therefore performed with periodic boundaries in 
z, assuming translational invar iance in x. The numerical 
counterpart of the fully developed velocity pattern shows 
good correspondence (Fig. [3h). 

Associated with the interfacial undulations are 
counter-rotatory streamwise vortices seen numerically in 
the secondary flow field Vy^i/z (Fig.[3]i). Experimentally, 
focal plane constraints make measurement of Vy, and so 
direct comparison with this numerical map impossible. 
However an analysis of Vz (not shown) indicates spatial 
modulation at an identical wavelength to that in Vx^ at 
least consistent with streamwise recirculation. 

We analyse the developed pattern by measuring the 
wavelength 2i\ j qz of the interfacial undulations across a 
range of G. We compare the size of the wavevector (f^ 



with linear stability analysis and full non-linear simula- 
tions of the final steady state in the DJS model (Fig. 
4). For the range explored we find a decreasing wavevec- 
tor, and good quantitative agreement between theory and 
experiment (higher values of G for which saturation of 
wavevector is predicted numerically are unattainable in 
our experiments). The amplitude A of the interface posi- 
tion undulations in y (Fig. 4, inset) increases with G and 
is a few percent of the channel height. 

Conclusion — Careful 3D velocimetric study of a shear 
banding wormlike micellar surfactant solution in mi- 
crochannel flow reveals a strongly enhanced confinement 
effect induced by the edges, and an interfacial instability 
with wavevector in the vorticity direction. Striking quan- 
titative agreement is obtained with a numerical study of 
the DJS model, without free parameter fitting. Given 
that this model is highly phenomenological, with only 
crude representation of the microscopic rheology, this 
suggests robustness of the phenomena to underlying mi- 
croscopic details, hinging only on the presence of a non- 
monotonic constitutive curve and a normal stress jump 
across the interface between the shear bands. The con- 
finement effect, here due to inhomogeneous banding on 
the walls, may prove generic in a broader class of inho- 
mogeneous confined flows, a conjecture we shall explore 
more fully in future work. 

Our observation of the interfacial instability in planar 
flow, combined with its linearly unstable nature, provides 
strong evidence in favour of scenario (ii) of the intro- 
duction. Experimental resolution unfortunately averages 
possible qx modes, predicted by linear stability analysis 
to be concurrent with qz ones (ongoing work). It thus 
remains an open challenge both theoretically and experi- 
mentally to delineate the effect of nonlinear coupling be- 
tween these modes; and to study the interplay between 
the interfacial instability and the enhanced confinement 
effect reported here. 
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